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Lattice QCD at finite density 

A primer 

Maria-Paola Lombardo 
Istituto Nazionale di Fisica Nucleare, Italy 

\ QCD at finite density presents specific challenges to lattice gauge theory. Nonetheless, 

a region of the QCD phase diagram up to moderately large baryon chemical potentials has 

1 been successfully explored on the lattice and new results and idea are continuously emerging. 

' I will outline the lattice formulation of QCD, introduce the calculational schemes cur- 

rently used to treat a nonzero baryon density, and mention lattice methods alternative to 
MonteCarlo, including the strong coupling expansion which might give access to the the 
I ' superconducting phase of QCD. The results for the critical line, and the different phases 

will be discussed highlighting the strength of the different methods, as well as the possible 
comparisons with phenomenological models. 

i— I , 

' 81. Introduction 

>: 

, Lattice discretization combined with importance sampling affords the possibility 

of doing first principles calculations of the properties of strongly interacting matter. ^ 
Questions to be addressed obviously include the study of QCD with physical values 
of the quark masses, and in the thermodynamic region within the range of current, 
and planned, experiments at BNL and CERN. 

While addressing these phenomenological points, one is lent to consider more 
theoretical questions - patterns of chiral symmetry, mechanisms of confinement, 
topological structures, gauge field dynamics...- on general grounds. To this end, it 
is not only legitimate, but also useful (and sometimes mandatory), to consider QCD 
in a larger parameter space. This can be achieved in a numerical simulation, where 
not only we can tune experimental parameters such as baryon chemical potential, 
^ ■ temperature, isospin chemical potential, but we can also play with the number of 

colors and flavors, the value of the bare masses, the gauge coupling, or we can make 
imaginary some external fields. In the following QCD will thus be a generic name for 
any of these variants - to be specified whenever needed - of the "real world" theory 
of strong interactions. 

In this note we will be mostly concerned with the phases of QCD at finite 
baryon density. Let me then mention, as a last introductory remark, that there are 
strong differences between physical mechanisms of phase transitions in QCD at high 
baryon density and high temperature. At high temperature the chiral transition is a 
transition from an ordered to a disordered state, characterized by light baryons; de- 
confinement is associated with string breaking due to recombination with light pairs; 
topological structures include instanton molecules. At high baryon density instabil- 
ities at the Fermi surface induce unusual patterns of chiral symmetry such as color 
superconductivity or superfluidity; string breaking might be further enhanced by 
recombination with real particles; topological structures are mostly likely instanton 
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chains. These differences between high temperature and high density might provide 
further insight into chiral symmetries, gluon dynamics and topological structures, 
and further motivate the study of a nonzero baryon density on a lattice. 

§2. Formulation 

Let us remind ourselves how to introduce a chemical potential fi for a conserved 
charge TV in the density matrix p in the Grand Canonical formalism, which is the 
one appropriate for a relativistic field theory: ^ 

p = e -{H-»N)/T (2 . 1} 

Z(T,fi) =Trp = J dcf)d^e- s ^ (2-2) 

The path integral representation of the grand partition function Z in the Eu- 
clidean space gives the temperature as the reciprocal of the imaginary time: 

i It 

S(<l>,il>) = J dt J d d x£{(j),ij) (2-3) 

with periodic boundary conditions in time for bosons (f>(t = 0, x) = <p{t = 1/T,x) 
and antiperiodic for fermions ip(t = 0, x) = —ip(t = 1/T,x). 

All in all, Z at finite temperature T and density fi is the partition function of a 
statistical system in d+1 dimension, where T is the reciprocal of the imaginary time, 
and \x couples to any conserved charge. This representation, which is the starting 
point for a lattice calculation, allows us to deal with thermodynamics and spectrum 
exactly on the same footing. 

The theory is regularised on a space time lattice: a regular four dimensional grid 
with N s points in each space directions, Nt points in the imaginary time direction, 
and spacing a. We refer to the very many excellent reviews and textbooks for 
background material on lattice field theory, and we briefly summarize here the specific 
aspects of lattice QCD thermodynamics which will be useful in the following. 

The temperature T on a lattice is the same as in the continuum: T = 1 /Nta, N^a 
being the lattice extent in the imaginary time direction (while, ideally, the lattice 
spatial size should be infinite). A lattice realisation of a finite density of baryons, 
instead, poses specific problems: the naive discretization of the continuum expression 
/z^7oV' would give an energy e oc ^ diverging in the continuum (a — ► 0) limit. 2 ) 

The problem could be cured by introducing appropriate counterterms, however 
the analogy between \x and an external field in the t h (temporal) direction offers 
a nicer solution by considering the appropriate lattice conserved current. 2 ) This 
amounts to the following modification of the fermionic part of the Lagrangian for 
the (^direction L° F : 

L° F (p) = Yw^+e " VWoToe"^ (2-4) 
while the remaining part of the Lagrangian is unchanged. This yields the current: 
Jo = -d^L = -d^L F °(p) = V^oe^+o + VW67oe- Ma ^ (2-5) 
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This representation of Jo is amenable to a simple interpretation: the time forward 
propagation is enhanced by e^ a , while the time backward propagation is discouraged 
by e~^ a ; hence, the link formulation generates a particles-antiparticles asymmetry. 
In addition, note that J* Jo = N — N as it should. An alternative way to look 
at the link formulation introduces an explicit dependence on the fugacity e^l T via 
an unitary transformation for the fields. 3 ) In this way L(p) = L(0), and the n 
dependence is on the boundaries, via the fugacity e^ T : iP(x + Nt) = —e fMaNT i{j(x) = 
—e fM ^ T ip(x). This is analogous to the continuum case. 4 ) 

§3. Calculational Schemes 

Having set up the formalism, the task is to compute 

Z = J dUd4,e- s{u ^ (3-1) 

where from now on the Lagrangian defining the Action will be that of lattice QCD, 
containing gluon fields U and quark fields ip. 

We have two options. We might integrate out gluons first: 

J dUdil>di>Z(T,n,4>,il>,U) ~ J di/>dij>Z(T,ii,ii>,il>) (3-2) 

This produce an effective approximate fermion model: the procedure is physically 
appealing, but not systematically improvable, but for one special (lattice) case (see 
below). Alternatively, we might integrate out fermions exactly, by taking advantage 
of the bilinearity of the fermionic part of the Lagrangian L = Lym + Lp = Lym + 

J dUdipd$Z(T,n,$,ip,U) = J dUe-( SYM W- lo9( - detM » (3-3) 

The "effective" model we build this way is exact: the price to pay being that its 
physical interpretation is not as clear as for effective fermion models. Anyway, this 
expression is the starting point for numerical calculations: the fact that in many 
cases they are highly successful tell us that the configuration space is well behaved 
enough that only a minor subset of configurations, although carefully chosen via 
importance sampling, suffice to produce reasonable results. 

3.1. Effective Fermionic Models: analytical approaches 

Let us start by following the first idea, namely integrating out the gluon fields 
so to define an effective fermionic Action. This is a time honored approach, leading, 
for instance, to the instanton model Hamiltonian, hence to the exciting discoveries 
on the QCD phase diagram of the last five years 5 ) . 

On the lattice, one very interesting approach leading to a fermionic model is 
provided by the strong coupling expansion: in the infinite gauge coupling limit the 
Yang Mills term decouples from the Action, and the integral over the gauge fields 
can be carried out exactly. 
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The starting point is the QCD lattice Lagrangian: 

3 

S = -1/2 E m( x )\x{x)U 3 (x) X {x + j) - x(x + j)U}(x) X (x)) (34) 

X j = l 

-1/2 ^ %(aO[x(s)E/o(s)x(s + 0) - x(x + 0)C/ t (x) X (x)] 

4 

-1/3 ^6/ 5 2 ^ [l-reTrt^(x)] 

x /J-,v=l 

+ ^2mxx 

X 

The X) X are the staggered fermion fields living on the lattice sites, the IPs are 
the SU(N C ) gauge connections on the links, the 77's are the lattice Kogut-Susskind 
counterparts of the Dirac matrices, and the chemical potential is introduced via 
the time link terms e M , e _M as discussed above. This time we have written down 
explicitly the lattice Action to show that the pure gauge term 

Sq = — 1/3 Ylx ^ = i4 [1 — reTrU^ u (x)] contains the gauge coupling in the 

denominator, hence it disappears in the infinite coupling limit. Consequently, one 
can perform independent spatial link integrations, leading to 

Z = j Yl dU t dxdxe- 1/4N ^<™>* {x)x{xmy)x{y) e- St (3-5) 

timelinks 

where ^2 <x y> means sum over nearest neighboring links, terms of higher order have 
been dropped, and we recognize a four fermion interaction. 6 ) Further manipulations 
yield the mean field effective potential: 

V eff (< ^ >, n) = 2cosh(rN t N c fi) + sinh[(N t + 1)N C < ^ >}/sinh(N t < ^ >) 

which we quote for further reference. A standard analysis of V e ff finally gives the 
condensate as a function of temperature and density, and allows the reconstruction 
of the phase diagram. 

More recently this approach has been furthered both in two 8 ) and three colors, 9 ) 
and new developments on cluster algorithms have appeared as well. 10 ) 

In order to describe in detail the rich physics of the finite density phase, one needs 
both to include higher order terms into the strong coupling expansion, as well as to 
go beyond a simple mean field analysis, which assumes an homogeneous background. 
The question is as to whether such improved strong coupling approaches would be 
able to generate a four fermion term with the correct flavor structure as well as order 
of magnitude, thus opening the possibility of a systematically improvable approach 
to finite density QCD, including the study of the superconducting phase. 

3.2. Effective Gluonic Models: Importance Sampling and the positivity issue 
Let us write again 

Z(T, (j) = [ dUe- {SYM{u ^- lo ^ dct M )) (3-6) 
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When det M > the functional integral can be evaluated with statistical meth- 
ods, sampling the configurations according to their importance (Sym(U) — log(det M)) 
For this to be possible the would-be-measure (det M) has to be positive. 

Let me mention at this point that the factorization method 7 ) might alleviate 
the problems of complex measures by guiding the simulations along a sensible path 
in the phase space. I will not dwell on this interesting development which is not 
really in the scope of an introductory review, but I wish to call on it the attention of 
the interested reader, as it really seems to offer some promise, and has been already 
tested in random matrix models. 

In QCD with an even number of flavors, and zero chemical potential, standard 
importance sampling simulations are possible if det M is real, which is true if = 
—PMP~ 1 where P is any non singular matrix. In the most popular lattice fermion 
formulation this holds: for Wilson fermions P = 75 and for staggered fermions 
P = I (note that this basically expresses a particle—antiparticle symmetry) . We will 
consider staggered fermions from now on. 

Consider now the relationship M\hb) = —M(—[ib) implying that reality is 
lost when Re^ / 0: the reality of the determinant is lost, and with it the possibility 
of doing simulations with non zero chemical potential, when we want to create a 
particle antiparticle asymmetry. On the other hand a purely imaginary chemical 
potential does not spoil the reality of the determinant: indeed, even if an imaginary 
chemical potential can be used to extract information at real chemical potential, it 
does not create any real particle-antiparticle asymmetry and it is natural that the 
fermion determinant remains real. 

Note that in QCD with two color the determinant remains positive with nonzero 
real chemical potential: indeed, in that case quarks and antiquarks transform under 
equivalent representation of the color group and are, essentially, the same particle. 
Other important models with a real determinant include finite density of isospin 12 ) 
and four fermion models. n ) These aspects will be reviewed by Don Sinclair at this 
meeting, 12 ^ so I shall not further discuss them. 

All in all, if we want to extract information useful for QCD at nonzero baryon 
density by use of standard MonteCarlo sampling we will have to use information 
from the accessible region: 

Re fi = 0, Im jjL < 
§4. Overview of the methods 

To begin with, it is useful to think of the theory in the T, plane. Let us 
then discuss the phase diagram from the perspective of analyticity and positivity of 
the partition function and of the determinant. One important consideration to keep 
in mind: the Gran Canonical partition function has to be positive. It is only the 
determinant which can change sign, or even be complex, on single configurations. 

Let us consider a mapping from complex fi to complex /i 2 . Because of the sym- 
metry properties of the theory, this mapping can be done without loss of generality. 
Let us note then that Z{n 2 ) is real valued for real /i 2 : this is a situation familiar 
from condensed matter: the partition function is real where the external parameter 



6 



M.-P. Lombardo 



is real, complex otherwise. 

The reality region for the partition function represents states which are physically 
accessible. The reality region for the determinant represents the region which is 
amenable to an importance sampling calculation: Kefx 2 < 0. The methods which 
have been applied so far are 

• fx = Derivatives, Reweighting, Expanded reweighting 

• fj 2 < Imaginary chemical potential 

4.1. Derivatives at fx = 0.0 

This is one early attempt at explor- 
ing the physics of nonzero quark den- 
sity: the derivatives can be formally 
computed at \x = 0. 13 ) The obvi- 
ous limitation is that we do not really 
know how far from the fx = axis can 
we get. Nonetheless, such derivatives 
are interesting per se, and the region 
where derivatives are clearly different 
from zero is the natural candidate for 
the application of other methods. 

4.2. Reweighting from fx = 

Back in the 80's Ian Barbour 
and collaborators proposed to calculate 
Z(jx) from simulations at fx = 0: 




Z 



I |M(/x)| 
\\M( f x = 0)\ 



(4-1) 



/i=0 



Fig. 1. Sketch of the phase diagram in the 
, T plane: the solid line is the chiral tran- 
sition, the dashed line is the Roberge Weiss 
transition. Simulations can be carried out 
at fi 2 < and results continued to the 
physical domain n > 0. The derivative 



and reweighting methods have been used 
so far to extract informations from simula- 
tions performed at p, = 0. The imaginary 
chemical potential approaches uses results 
on the left hand half plane. Different meth- 
ods could be combined to improve the over- 
all performance. 



In other words, the chemical potential fx 
of the target ensemble at that of the sim- 
ulation ensemble - fx = - are different: 
the properties of the target ensemble 
can be inferred from those of the sim- 
ulation ensemble, provided that there is 
a sizable overlap between the two. 14 ) 

At T = the Glasgow procedure 
fails because of a poor overlap (aside, the strong coupling calculations were quite 
useful to asses these problems) , and it is instructive to study the overlap problem as 
seen in the Gross Neveu model, where there is no sign problem, u ) and the results 
obtained with reweighting methods can be compared with those of exact simula- 
tions. 15 ) 

The distribution of the order parameter (the a particle) helps visualizing the 
problem (see Fig. 2): the order parameter distributions in the two phases do not 
overlap. 11 ^ 

It is interesting to note that there are indeed examples of successful reweighting 
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<c> for (l upd =0.0 



0.40 0.45 0.50 0.55 0.60 0.65 0.70 



0.75 0.80 0.85 0.90 



<0> for H„pd= - 625 



~l~rh-m 



0.10 0.20 0.30 0.40 0.50 0.60 



Fig. 2. The overlap problem in the Gross Neveu model: 5 ' the exact distributions of the < a > 
fields in the broken (left) and symmetric (right) phase are non-overlapping. 



at fx 7^ 0. Let us consider 1-dim SU(3) which can be exactly solved at nonzero 
baryon density 16 ) and can thus serve as a test bed for experiments. The distribution 
of the partition function zeros was computed and found to reproduce the correct 
results one the statistics were sufficiently high. 22 ) 




Fig. 3. Successful reweighting in one dimensional QCD: ' exact Z's zeros in the complex plane 
(diamonds) and the cloud of zeros obtained from reweighting with a very poor statistics(leftt). 
On the rigth the zeros from an high statistics reweighting simulation reproducing the exact 
result. 

The conclusion from these early studies was that reweighting fails in QCD at 
zero temperature because of a poor overlap, and that the reason behind the failure is 
practical rather than conceptual: the situation can be ameliorated if a better starting 
point were used. 

4.3. Fodor and Katz's multiparameter reweighting 

The prescription for ameliorating the overlap is due to Fodor and Katz 17 ) 18 ) 
whose Multiparameter reweighting use fluctuations around T c at \i = to explore 
the critical region. Making reference to Fig. 2, and oversimplifying: instead of 
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trying to reweight the distribution at zero temperature in the broken phase, which is 
obviously hopeless, one might hope that a distribution generated at zero density, and 
close to the critical temperature, bears more resemblance with the target distribution 
along the critical line, and is thus amenable to a successful reweighting. 

The strategy was applied to QCD 17 ) 18 ) . The improvement obtained is impressive 
and produced the first quantitative results for the critical line at nonzero chemical 
potential in QCD: we will come back to this in the section on results. A multistep 
reweighting proposed by Crompton 19 ) might well produce a further improvement. 

4.4. Taylor Expanded Reweighting 

The Bielefeld-Swansea collaboration suggested a Taylor expansion of the reweight- 
ing factor as a power series in A = li/T, and similarly for any operator 21 ) 2 °) . 

This strategy is computationally very convenient as it greatly simplifies the 
calculation of the determinant. Expectation values are then given by 

, n , _ ((Oo + giA + O2X 2 + ■■■) exp(7^A + Tl 2 \ 2 + ...- AS g )) x=Q>Po 

{ kM (exp(K 1 \ + K 2 \* + ...-AS g )) x=0 , (h • [ ' 

Results - to be discussed later- have been obtained both for the critical line and 
thermodynamics. 

4.5. Imaginary Baryon Chemical Potential 

This method uses information from all of the negative 11 2 half plane (Fig. 1) to 
explore the positive, physical relevant region. An imaginary chemical potential v in 
a sense bridges Canonical and Grand Canonical ensemble: 23 ) 

Zc{M) = 7^j Q dvZ Q c{iv)e-^ N (4-3) 

The main physical idea behind any practical application is that at fx = fluctuations 
allow the exploration of iV& 7^ hence tell us about /z ^ 0. Mutatis mutandis, this 
is the same condition for the reweighting methods to be effective: the physics of the 
simulation ensemble has to overlap with that of the target ensemble. 

A practical way to use the results obtained at negative fx 2 relies on their ana- 
lytical continuation in the real plane. For this to be effective 22 ) T) must be 
analytical, nontrivial, and fulfilling this rule of thumb: 

X (T, ll) = 8p(n, T)/dii = d 2 logZ(ii, T)/ V > (4-4) 

This approach has been tested in the strong coupling limit 22 ) of QCD, in the 
dimensionally reduced model of high temperature QCD 24 ) and, more recently, in the 
two color model. 25 ) 

Results (to be discussed later) have been obtained for two, 26 ) four 27 ) and three 28 ) 
staggered flavors. 

§5. Results 



The methods just outlined above are workarounds, not real solutions: they are 
practical tools to circumvent a problem, and, as such, it it not surprising that they 
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have to be applied with a grain of salt, and that their performance depends on the 
thermodynamic region which is being explored. 

In this section I will go through the main issues which have been addressed so far: 
the critical line, the hadronic phase, the "Roberge Weiss" regime, the quark gluon 
plasma phase, highlighting the main strengths of the various methods alongside with 
the results. 

5.1. The Critical Line 

The critical line has been obtained either by Fodor and Katz 17 ) 18 ) and by the 
Bielefeld Swansea collaboration within the multiparameter reweighting or the ex- 
panded reweighting approach, which gives Z{y) from Z(0). The location of the end 
point follows naturally within this framework, and its first determination was given 
in 18 ). 

De Forcrand and Philipsen have also noticed that the analytic continuation of 
the critical line from an imaginary /i is possible, 26 ) and have indicated and applied 
a strategy for the location of the end point. 28 ) 

Also for the calculation of the critical line the consideration of the T, fi 2 plane 
helps the analysis. 27 ) Model analysis suggests the following parametrization, con- 
firmed by numerical results: 



(T + aT c )(T - T c ) + kfi 2 = 0, k > 



(5-1) 



It encodes reality for real fi 2 , contains the physical scale T c , is dimensionally con- 
sistent, gives T{n = 0) = T c , T(fi ^ 0) < T c . For instance, the second or- 
der approximation T(T — T c ) + /u 2 /(8/n2) = to the Gross Neveu Model criti- 
cal line: 1 — fi/ Uq = 2T/Soln(l + e _/i//T ) is good up to \i ~ T c , and from the 
plot in Fig. 4 can see that a second order approximation is good up to \i ~ T c 



Fig. 



GN critical line 
0(mu ) approximation 



Studies of the critical line have 
indeed found that a simple polyno- 
mial approximation suffices to describe 
the data, within the current precision. 
Progress on the precision is demon- 
strated in Fig. 5 28 ) where the new re- 
sults on the three flavor model 28 ) are 
superimposed to the older ones with 
two 26 ) and four flavors. 27 ) 

A crucial issue remains the deter- 
mination of the endpoint, for which 
the first estimate was given within 
the reweithgting method Te = 160 ± 
3.5MeV, n E = 725 ± 35iWW. 18 ) Re- 
sults with improved precisions show a dependence on the mass values. This makes 
mandatory an extrapolation to physical values of the quark masses, which , in turn, 
implies a good control on the continuum limit. 



0.6 
mii/SO 



4. The critical line of the 3d Gross Neveu 
Model and its polynomial approximation. 
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5.2. Hadronic Phase T <T r I: the hadron resonance gas model 

In this region observables are a con- 
tinuous and periodic function of fii/T , 
analytic continuation in the [i 2 > half 
plane is always possible, but interesting 
only when x<?(m = ®>T) > 0. 

Taylor expansion and Fourier de- 
composition are natural parametriza- 
tion for the observables. 27 ^ In partic- 
ular, the analysis of the phase diagram 
in the temperature-imaginary chemical 
potential plane suggests to use Fourier 
analysis for T < T c , as observables are 
periodic and continuous there. Note 
that in the infinite strong coupling all 
of the Fourier coefficients but the first 
ones will be zero (cfr eq. 3-6). 
As the chiral condensate is an even function of the chemical potential, its Fourier 
decomposition reads: 

{W) = J2 a F cos(nN t N cf Mi) (5-2) 

n 

which is easily continued to real chemical potential. 

One cosine fit is actually enough to describe the data up to T ~ T c in the four 
flavor model: 27 ' adding a term cas(2NtN c /j,) in the expansion does not modify the 
value of the first coefficients and does not particularly improve the x 2 ■ 

This phase was also studied within the reweighting approach. 29 ) It was confirmed 
that one single hyperbolic cosine is an excellent approximation to the data, and the 
result has been interpreted within the framework of the hadron resonance gas model, 
whose partition function has the single hyperbolic cosine form as the one given by 
the strong coupling expansions eq. 3-6. 

5.3. The Hadronic Phase II: The order of the phase transition, and related endpoints 

The analytic continuation of an observable O is valid till \i < /i c (T), where the 
critical value (J* C (T) has to be measured independently. The value of the analytic 
continuation 0(fj, c ) of an observable O at [x c defines the discontinuity at the critical 
point. In turns, this allows the identification of the order of the phase transition. 
One might wonder which is the meaning of the analytic continuation for (i > fj, c (T), 
the one which we have to chop by hand. 

It is natural to interpret such analytic continuation as the metastable branch of 
the observable we are considering, for instance < tptp > : it follows the secondary 
minimum of the associate Landau Ginsburg potential and determines the spinodal 
point \i* according to < -tptp >= A(fi—fi*)P. The discontinuity < -iptp > (// c ) is related 
to (n c — //*), and both shrinks to zero at the endpoint of a first order transition. 

The analytic continuation of the results in the hadronic phase, when cross exam- 




0.9 - 



0.5 1 1.5 2 2.5 3 

~B c 

Fig. 5. Summary plot 28 ' for the critical line 
for N f = 2, 26) N f = 3, 28) N f = 4 27) from 
imaginary chemical potential calculations. 
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ined with the results for the critical line, offers an alternative way to study the order 
of the phase transition: the transition will be second order is the zero of the analytic 
continuation matches the critical point, first order otherwise. The endpoint of a first 
order transition can be detected by monitoring (/i — /x*) as a function of temperature. 

Clearly this approach is specific to 
the imaginary chemical potential calcu- 
lations, and give results which should be 
cross checked with others. For the time 
being it was confirmed 27 ^ that the tran- 
sition in the four flavor model remains 
of first order at nonzero density. 

5.4. The Roberge Weiss Regime: Tq < 
T<T E 

Let us consider region which is com- 
prised between the deconfinement tran- 
sition, and the endpoint of the Roberge 
Weiss transition: the analytic continua- 
tion is valid till fi = oo but the interval 
accessible to the simulations for fi 2 < 
is small, as simulations in this area hits 
the chiral critical line. 

The bright side of this is that the 
nature of the critical line can then be 
studied without need for an analytic 
continuation. In Fig. 3 we show the 
clear correlation between the Polyakov 
Loop and the chiral condensate at fi = 
0.15 27 ). 33 ) 

The correlation between chiral and deconfining transition persists at nonzero 
imaginary chemical potential, see Fig. 727) 33) Obviously, these observations can 
be immediately continued to real chemical potential : if the difference between the 
value of the critical for confinement and chiral symmetry is zero in a finite interval 
within the analyticity domain (in our case, for fi 2 < 0), it will be zero everywhere 
within the same domain. 

It is also of interest to note that the non-applicability of perturbation theory in 
this region is almost a theorem: indeed the analytic continuation of the polynomial 
predicted by perturbation theory for positive fi 2 would never reproduce the correct 
critical behavior at the second order phase transition for /j, 2 < , and it is then ruled 
out. 

5.5. The QGP phase : T E <T 

At high temperature, in the weak coupling regime, perturbation theory might 
serve as a guidance, suggesting that the first few terms of the Taylor expansion 




Fig. 6. The Bielefeld-Swansea collaboration 
data contrasted with the hadron resonance 
gas model. 29 - ) 
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might be adequate in a wider range of chemical potentials. This confirms that 
the Roberge Weiss critical line has to be strongly first order at high temperature. 



Several analytic models have been 
proposed to describe the properties of 
this phase. As one example, we show 
how a QGP liquid model prediction de- 
veloped by Letessier and Rafelsky: 31 ) 

AP/T 4 = 2(1 -2a s /^)(WT) 2 ) (5-3) 



, 1 i 



where 



com- 



f 



Fig. 7. Correlation between < ipip > and 
Polyakov loop at m = 0.15. 33) The ob- 
served correlation of chiral and deconfining 
transition must persist at nonzero baryon 
density. 27 - 1 



Us = a s {2^J (7rT) 5 , rqj 

pares favorably with reweighting data 
by Z.Fodor, S.Katz and K.K.Szabo. 30 ) 
Very similar conclusions were reached 
by analyzing quasi-particle models. 32 ) 

As a further practical tool for the 
analysis of this phase, it was pro- 
posed 33 ) to define an "effective prefactor 
plot" APK L(Nt=4) /(T 4 2(fi/T) 2 ) which 
should equal, for instance 2(1 — 2a s /ir) in the model discussed above, and, in general, 
which would serve to assess by eye the \i dependence, if any, of the prefactor of the 
quadratic term. 

It was found that the results approach the perturbative limit at large chemical 
potential, but corrections at small chemical potential are clearly visible, and become 
much less severe while approaching the large \i limit. It would be interesting to 
understand the connection between these numerical observations and the theoretical 
work on effective positivity on dense matter. 34 ) 



§6. Summary /Out look 

QCD at nonzero baryon density is now an active field of research with results 
emerging from different methods. 

The critical line for two, three, two plus one, and four flavors has been computed 
by various methods, with a substantial agreement. The critical line is well described 
by a polynomial, and this result can be interpreted in terms of simple models. 

The endpoint has been located by different methods, and its determination is 
currently being sharpened. In the four flavor model, when the transition is of first 
order, the chiral and deconfining transition remain correlated at nonzero chemical 
potential. 

The analytic continuation from an imaginary chemical potential gives access to 
the often evasive physics of the metastable branch. It might afford an alternative 
way to locate the end point and tricritical point. 

Three different regimes have been considered and discussed: the hadronic phase 
results are consistent with the hadron resonance gas model, both from reweighting 
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calculations and imaginary chemical potential approach; the Roberge Weiss regime 
is eminently nonperturbative, and in this regime we have the possibility to study 
the nature of the chiral transition at nonzero chemical potential without performing 
any analytic continuation; the Quark Gluon Plasma phase seems well described by 
simple models, but corrections are visible and still need being quantified. 

Assessing the validity of simple models at high temperature/density, besides 
being extremely interesting per se, would also open the possibility of doing non 
equilibrium calculations based on these models. The imaginary chemical potential 
approach seems to be particularly well suited for this task. 

The three methods we have dis- 
cussed are mature for quantitative stud- 
ies in realistic cases. A nice possibility is 
offered by a combination of these meth- 
ods: for instance either reweighting or 
direct calculations of derivatives could 
be performed at nonzero fi to improve 
the accuracy of the results at negative 
/x 2 , and of the ensuing analytic continu- 
ation to real [i. 

Finally, it might well be that other 
methods such as xQCD which in the 
past encountered difficulties at zero 
temperature and nonzero baryon den- 
sity will prove successful in the richer 
high temperature regime. 12 > 




T/T„ 

Fig. 8. Lattice results by Fodor, Katz and Sz- 
abo contrasted with a Quark Gluon Plasma 
liquid model. 31 ' 
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